Measurement of the CKM angle a {(^2) with B ^ pp decays at Belle and BaBAR * 
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We overview recent measurements m B ^ pp decays which are based on data samples collected 
at the PEP-II and KEKB asymmetric-energy e^e~ colliders with the BABAR and Belle detectors. 
Special emphasis is given to the determination of the CP-violating coefficients A and S from an 
analysis of B'^ —* p^ p~ decays. The values of A and S, branching fractions, and longitudinal 
polarization fractions of B pp decays are used to constrain the Cabibbo-Kobayashi-Maskawa 
phase a{(f)2) using an isospin analysis; the solution consistent with the standard model is 71° < 
q(<?!>2)<113° at 68% C.L. 



o 
o 



00 



> 

o 
cn 
o 
l> 
o 

Oh: 



X 



I. INTRODUCTION 

CP violation in the Standard Model can be explained 
by the presence of an irreducible complex phase in the 
Cabibbo-Kobayashi-Maskawa [1] (CKM) quark-mixing 
matrix. The unitarity of the CKM matrix leads to six 
triangles in the complex plane. One such triangle is 
given by the following relation among the matrix ele- 
ments: V*f,Vud + V*i,Vcd + V^lVtd = 0. The phase angle 
(j)2^, defined as arg[-(VtdV*i^)/{VudV can be deter- 
mined by measuring a time-dependent CP asymmetry in 
b — *■ uud decays such as — j-tt+tt^, p^n^, and p'^ p~ ■ 
The time-dependent decay rate for B P~ decays 

tagged with B^{q = 1) and B° {q = —1) mesons is given 
by 



.-|At|/r^o 



4r 



-{1 + q[App cos(AmAt) (1) 
+Spp sin(ATOAi)]}, 



where tbo is the B^ lifetime, Am is the mass difference 
between the two B'^ mass eigenstates. At — tcp — ttag, 
and App {BaBAR's definition is Cpp = —App) and Spp are 
CP asymmetry coefficients to be obtained from a fit to 
the experimental data. If the decay amplitude is a pure 
CP -even state and is dominated by a tree diagram, Spp — 
sin(2(^2) and App = 0. The presence of an amplitude 
with a different weak phase (such as from a "penguin" 
diagram) gives rise to direct CP violation and shifts Spp 
from sin(2^2)- However, the size of the loop amplitude is 
constrained by the branching fraction of B^ p'^p'^ [2], 
indicating that this effect is small. 

The CP- violating parameters receive contributions 
from a longitudinally polarized state (CP-even) and 
two transversely polarized states (an admixture of 
CP-even and CP-odd states). Recent measurements 
of the polarization fraction by Belle and BaBAR [3, 4] 
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show that the longitudinal polarization fraction is 
approximately 100% (/l = 0.968 ± 0.023 [5]). /l can 
be extracted from a fit to helicity-angle distribution. 
The angular decay rate (Pr/{T dcos9+ dcos9-) is given 
by I {/l cos^ 0+ cos2 0_ + i(l - h) sin^ 9+ sin^ 9. } , 
where, 9^ is the angle between the direction of the tt^ 
from the p^ and the negative of the B^ momentum in 
the p^ rest frame. 

II. MEASUREMENTS 

The common features oi B ^ pp analyses are: (1) 
relatively small signal yields; the branching fractions 
of S — > decays are in the order of 10~^ — 10~^, 
(2) large width of p mesons results in the large back- 
ground; the fraction of signal events in most analysis is 
less than 1%, (3) there are several background sources: 
e'^e~ qq {q = u,d,s,c) continuum, 6— > c, and 6— > u 
backgrounds, (4) significant amount of events with mul- 
tiple reconstructed B candidates, (5) various variables 
are used in the likelihood functions to distinguish among 
signal and backgrounds. 

Both Belle and BaBAR analyses identify B ^ pp de- 
cays using the beam-energy constrained mass M^^ = 

V -^bcam ~ P% (called as beam-energy-substituted mass, 
TTiES: in BaBAR) and energy difference AE = Eg — E^^^^, 
where -Ebcam ^be beam energy, and and p\ are the 
energy and momentum of the reconstructed B candidate, 
all evaluated in the center-of-mass (CM) frame. 

The dominant background originates from e+e^ — > 
99 (9 = u,d,s,c) continuum events. To separate qq 
jet-like events from spherical-like BB events. Belle uses 
event-shape variables, specifically, modified Fox- Wolfram 
moments combined into a Fisher discriminant [6], and 
9g, the polar angle in the CM frame between the B di- 
rection and the beam axis. The Fisher discriminant and 
9g are used to form a ratio of signal and background 
likelihoods TZ. In the B^ P^P^ analysis Belle also 
requires |cos0t| < 0.8, where 9t is the angle between 
the thrust axis of the candidate tracks and that of the 
remaining tracks in the event. In the BaBAR analyses qq 
background is suppressed by requiring | cos^tI < 0.8 and 
making use of a neural network discriminant J\f which is 
based on several event-shape variables. Bellow we de- 
scribe B ^ pp measurements in detail. 
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TABLE I: Reconstruction requirements used in 
analysis. The barrel and (endcap) regions in the electromag- 
netic calorimeter at Belle are defined as 32° <6 < 129° and 
(17° < 61 < 32° and 129° < 6* < 150°), where 8 denotes the polar 
angle with respect to the beam axis. 



Cut 


BABAR 


Belle 


(MeV) 


50 


50 (barrel) 90 (endcap) 


M-,^ (MeV/c") 


100. - 160. 


117.8 - 150.2 


M,±,o (GeV/c^) 


0.5 - 1.0 


0.62 - 0.92 


cos6l-j- 


-0.9 - 0.98 


-0.8 - 0.98 


Mt, (GeV/c^) 


5.25 - 5.29 


5.23 - 5.29 






5.27 - 5.29 {CP fit) 


AE (GeV) 


-0.12 - 0.15 


-0.2 - 0.3 






-0.12 - 0.08 (CP fit) 



^ BaBar 









5.25 5.26 5.27 5.28 5.29 
iTigs (GeV/c-) 
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FIG. 2: BABAR the At distributions of a signal-enriched sam- 
ple for (a) B'^ and (b) B tagged events, (c) shows raw asym- 
metry. The dashed lines denote qq background, the dotted 
lines are the sum of background, and the solid lines are the 
total. 
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FIG. 1: BABAR. (a) tubs, (b) AE, (c) cosO, and (d) m^±^o 
for the highest purity tagged events. The curves show fit 
projections: the dashed line is background and the solid line 
is the total. 



B — > decays are reconstructed by combining 

two oppositely charged pion tracks with two neutral pi- 
ous. The mesons are selected combining tt^ with tt*^ 
candidates. The tt" candidates are reconstructed from 
77 pairs. Main event reconstruction requirements are 
listed in Table I. A flavor of the B meson accompany- 
ing the B^ P^ P^ candidate is identified via its de- 
cay products. Tagging algorithms yield the flavor of the 
tagged meson and a flavor-tagging quality. In events 
with multiple reconstructed B candidates the best can- 
didate is selected based on the tt'^ masses, i.e. mini- 
mizing X^ttJ (™77 ~ JTijro)^. The fraction of signal de- 
cays which have at least one tt^ track incorrectly iden- 
tified but pass all selection criteria is 13.8% and 6.5% 
for BaBAR and Belle, respectively. These are reffered to 
as "self-cross-feed" (SCF) events. The following compo- 
nents are distinguished in the analyses: signal and ptttt 
non-resonant decays, SCF events, continuum background 
[qq) , charm B background (& — > c) , and charmless (6 u) 
background. The b u background is dominated by 
B^{pn, a^TT, ai p, p^ p^) decays. 

BaBAR obtains the signal yield, /i, and CP-violating 



FIG. 3: Belle: projections in Mbc (left) and AE (right) for 
the high-purity tagged events of a sample enriched in signal. 
The curves show fit projections: dashed is p'^p'+pmr, dotted 
is qq, dot-dashed is b—>c, long-dashed is h^u, and solid is 
the total. 

parameters .Ap+p- andiSp+p- from an unbinned extended 
maximum-likelihood (ML) fit to mESi ^E, /At, m^±^o, 
cos 9± , and N distribution of 33902 events [7] . 

The Belle analysis is organized in two main steps: (a) 
we first determine the yields of signal and background 
components from an unbinned extended ML fit to the 
three-dimensional M^,^ - AE - TZ distribution, (b) we 
perform a fit to the At distribution of 18004 events to 
determine the CP parameters .4p+p- and Sp+p- The fit 
results are presented in Figures 1, 2, 3, and 4 and are 
listed in Table 111. The fraction of longitudinal polar- 
ization and fraction of non-resonant events (6.3 ± 6.7%) 
were measured in our previous analysis [3]. The branch- 
ing fraction, f^, and CP asymmetries measured by Belle 
are similar to those obtained by BABAR. The values 
.4p+p- and Sp+p- are also consistent with no CP vio- 
lation (^ = 5 = 0). 

B. B± ^ p±p° 

The main reconstruction features of the analysis are 
the same as those for the B° P^P - Event selection 
requirements are listed in Table II. 
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FIG. 4: Belle: the At distribution and projections of the fit 
for high-purity tagged events: (a) tags, (b) B tags. The 
raw CP asymmetry is shown in (c). The hatched region shows 
signal events. 
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FIG. 5: BABAR. (a) tues, (b) AE, (c) cos0, and (d) m^±^o 
distributions for a signal-enriched samples along with fit pro- 
jections: the dashed lines show qq and BB backgrounds, the 
solid line is the total. 



BaBAR obtained the yields of 



and 



B P fo decays, polarization, and charge asymme- 
try Ap±pO = {Nb- - Nb+)/{Nb- + Nb+) using an un- 
binned extended ML fit to tties, Ai?, At, m^±^o, to^+^- , 
cos^i, cos 6*0, and JV distribution of 74293 events [8]. 
The charmless b ^ u background is dominated by rj' p^, 
K*'^p^, fliTT^, afn^, afp'^, and a^p^ decays. 
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FIG. 6: Belle: AE (left) and _Mbc (right) distributions. The 
dashed curve is the sum of BB and continuum backgrounds, 
the dot-dashed curve is the signal, the hatched region shows 
BB, and the solid curve is the total. 



TABLE II: Reconstruction requirements used in B^ — > P^p'^ 
analysis. 



Cut 



BABAR. 



Belle 



E^ (MeV) 



50 



(MeV/c") 
M^i^o (GeV/c^) 

„2\ 



(GeV/c 



cos 9± 
cos 6{) 

(GeV/c^ 
AE (GeV) 



100. - 
0.396 - 
0.520 - 

-0.9 - 
-0.95 



160. 
1.146 
1.146 
0.95 
0.95 



'<^*^ (GeV/c) 



5.26 < 
-0.15 - 0.15 



50 (barrel); 100 (endcap) 
118. - 150. 
0.65 - 0.89 
0.65 - 0.89 



> 0.5 
5.272 < 
-0.4 - 0.4 





FIG. 7: BABAR niEs and AE distributions and projections of 
the fit for a signal-enriched sample. The dashed curve is full 
background, the small solid curve is signal, the dotted curve 
is B^p^ fo, and the large solid is the total. 

The signal yield in Belle analysis [9] is obtained from 
a fit to AE distribution. To measure the polarization, 
Belle bins in cos6'± and cosSq and determine the signal 
yield for each bin from the fit to AE distribution. The 
polarization is obtained from a simultaneous fit to two 
background-subtracted helicity-angle distributions. The 
results of the fits are shown in Figures 5 and 6 and listed 
in Table III. 



III. B° 



p'^p*^ decays and mea- 



BaBAR finds evidence for B'^ 
sures its branching fraction and polarization using a sam- 
ple of about 348 million BB pairs [2]. Events are se- 
lected from the region 5.24 GeV/c^ < Mj^^ < 5.29 GeV/c^, 
\AE\ < 85 MeV, and are required to satisfy 0.55 < 
M^+„^ < 1.05 GeV/c2, and |cos6io| < 0.98. In events 
with multiple B candidates one is selected based on the 
the best of a four-pion vertex fit. Additional suppres- 
sion of the dominant continuum background is achieved 
using the flavor tagging information. The data sample is 
divided into seven tag-quality intervals, ctag- The B^ — > 
pOpO event yield and polarization are obtained from 
an unbinned extended ML fit to to_bsi ^E, m^+^-j^ 2j 
cos6'i_2, A/", and Ctag distribution of 65180 events. The fit 
also allows to obtain the yields for B'^ p^ /o(990) and 
SO ^/o (980) /o(980) decays. The charmless background 
is dominated by B*^ — > o^-k^ events which number is a 
free parameter in the fit. Other b ^ u background modes 
include: B ^ {p° , p+ p^ , pir), and B° ^ p+ p^ . The 
fit results are shown in Fig. 7 and listed in Table III. 
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TABLE III: Summary of B—>pp measurements. 



Decay 


Quantity 


BABAR 




Belle 




HFAG 






Value 


Ns^g 


£(fb-i) 


Value 


Ns^g 


£ (fb-i) 


Value 




Br (10~^) 


23.5 ±2.2 ±4.1 


615 ±57 


316 


22.8±3.8+H 


194 ± 32 


253 


no 1 +3.2 
^■^•^-3.3 


B ^ p p+ 


U 


0.977 ± 0.024t^ j^la 


615 ±57 


316 


o.94it;;:oi^o ± 0.030 


194 ± 32 


253 


0.968 ±0.23 




0.07 ±0.15 ±0.06 


615 ±57 


316 


0.16 ±0.21 ±0.07 


372 ± 43 


492 


0.11 ±0.13 




Sp+p- 


-0.19 ±0.21^5;;°^ 


615 ±57 


316 


0.19 ±0.30 ±0.07 


372 ± 43 


492 


-0.06 ±0.18 




Br (10"'^) 


16.8 ±2.2 ±2.3 


390 ± 49 


210.5 


31.7±7.lt|:? 


58.7 ± 13.2 


78 


18.2 ±3.0 


B-^P^p' 


U 


0.905 ± 0.042j:^;[]2? 


390 ± 49 


210.5 


0.948 ±0.106 ±0.021 


58.7 ± 13.2 


78 


n 010+0 044 

U. JiZ_Q Q4g 




Ap±pO 


-0.12 ±0.13 ±0.10 


390 ± 49 


210.5 


0.00 ±0.22 ±0.03 


58.7 ± 13.2 


78 


-0.08 ±0.13 


B ^ p»pO 


Br (lO-*^) 
U 


1.16to.3B ±0.27 
0.86lE;;i:^ ± 0.05 


981^3? 


316 
316 




1.16 ±0.46 




30 60 90 120 150 180 

^2 (degrees) 

FIG. 8: 1 - C.L. vs. (j>2- The horizontal lines denote C.L. = 
68.3% (solid) and C.L. = 90% (dashed) 

IV. CONSTRAINT ON a (02) 

We constrain (/)2 using an isospin analysis [10], which 
allows one to relate six observables to six underlying pa- 
rameters: five decay amplitudes for B ^ pp and the an- 
gle 02- The observables are the branching fractions for 
B p+p-, p+p°, and pV" (listed in Table III) [5]; the 
CP parameters Ap+p- and Sp+p-; and the parameter 



ApOpO for B p^ p'^ decays. The branching fractions must 
be muhiphed by the corresponding longitudinal polariza- 
tion fractions (taken from Table III) [5] . We neglect pos- 
sible contributions from electroweak penguins and 1=1 
amplitudes [11] to B^ — » p'^ p~ ■ We follow the statis- 
tical method of Ref. [12] and construct a x^(</'2) using 
the measured values and obtain a minimum (denoted 
Xmin); then scan 02 from 0° to 180°, calculating the 
difference Ax^ = X^(02) ~ Xmin- We insert Ax^ into 
the cumulative distribution function for the x^ distribu- 
tion for one degree of freedom to obtain a confidence 
level (C.L.) for each 02 value. The resulting function 
1 — C.L. (Fig. 8) has more than one peak due to ambi- 
guities that arise when solving for 02- Because ApOpO is 
not yet measured, we allow this observable to float; this 
produces the "flat-top" regions in Fig. 8. The solution 
consistent with the Standard Model is 71° < 02 < 113° 
at 68% C.L. or 67° < 02 < 116° at 90% C.L. Recently, 
a different model-dependent approach to extract 02 us- 
ing flavor 5[/(3) symmetry has been proposed [13]. This 
method would give more stringent constraints on 02. 

In summary, we present recent measurements in _B — > 
pp decays. These measurements are used to constrain 
the angle 02 using an isospin analysis. 
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